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ABSTRACT 

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) 
observed eleven transits of TrES-3b over two years in order to constrain system 
parameters and look for transit timing and depth variations. We describe an 
updated analysis protocol for APOSTLE data, including the reduction pipeline, 
transit model and Markov Chain Monte Carlo analyzer. Our estimates of the sys- 
P_i ■ tem parameters for TrES-3b are consistent with previous estimates to within the 

2a confidence level. We improved the errors (by 10-30%) on system parameters 
^ j like the orbital inclination (iorb), impact parameter (6) and stellar density (p^,) 

compared to previous measurements. The near-grazing nature of the system, 
and incomplete sampling of some transits, limited our ability to place reliable 
uncertainties on individual transit depths and hence we do not report strong ev- 
idence for variability. Our analysis of the transit timing data show no evidence 
for transit timing variations and our timing measurements are able to rule out 
Super-Earth and Gas Giant companions in low order mean motion resonance 
with TrES-3b. 

Subject headings: eclipses, stars: planetary systems, planets and sateUites: fun- 
damental parameters, individual: TrES-3b 



1. Introduction 

When an extrasolar planet eclipses its host star, the event is referred to as a transit or 
primary eclipse. During a transit, observers can detect dips in starlight caused by the planet 
obscuring a portion of the stellar disk when it passes in front of the star. The transit method 
applies to those systems where the orbital inclination of a planet is close to 90° {i.e. edge- 
on) with respect to the o bserver's sky-plane. The first transit observations were made by 



Charbonneau et al.l ( 120001 ). As of August 2012, more than 200 planets (exoplanet.eu) have 
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been detected using the transit method. The search for new exoplanetary systems via transits 
has advanced to space- based missions yith the launch of the European S pace Agency's 



(ESA ) CoRoT sateUite (jPridlund et al.l l2006l ) and NASA's Kepler mission ( iBorucki et aL 



20101). The objective of transit detection and follow-up is mainly to catalog and improve 
measurements of system parameters. Studying the characteristics of extrasolar planetary 
systems is crucial for developing theories of planet formation that can adequately explain 
the origin and evolution of all planetary systems (including our own). 

The target discussed in this pap er, TrES-3b is a Hot- J upiter with one of the shortest 
orbital periods known (P=1.3 days lO' Donovan et al.l 120071 ) among the exoplanets. The 
planet orbits a G-type star (Tpf j = 5720K), and has a mass and radius of Mp =1.92 Mjup 



and Rp =1.29 -Rjup, respectively (lO 'Donovan et al.ll2007l ). Due to its large impact parameter 



TrES-3b's transit has more of a v-shape than the typical u-shape. However, the transit is not 
grazing and the ability of various observers to consistently measure the transit parameters 
seem to indicate that the planet's disk cornpletely enters the s tellar disk during the transit 
( lO'Donovan et al.ll2007l : ISozzetti et al.ll2009t iGibson et al.ll2009l ) . A blended eclipsing binary, 
which could account fo r the v-shape, has been ruled out from the radial velocity analysis 
( lO'Donovan et al.ll2007l ). However, the v-shape of the transit makes measurements of transit 
properties challenging at shorter wavelengths, where stellar limb-darkening further degrades 
the trapezoidal u-shape of the transit. The level of insolation received on the dayside of 
TrES-3b due to its proximity to its host star should place it in the class of warm Hot 



Jupiters with a temperatu re inversion in its atmosphere (iHubeny et al.ll2003l : iBurrows et al. 



20071 : iFortney et al.ll2008[ ). However, secondary eclipse observations fr om the ground and 
space have shown that t he evidence for an inversion layer is not strong (jde Mooij fc Snellen 
20091 : iFressin et al.ll2010l ). TrES-3b's lack of an inversion layer is consistent with the idea that 
UV radiation from chromospherically active stars (like TrES-3b's parent star) systematically 
destroys absorbers in the upper planeta ry atmosphere, thereby preventing the formation of 
an inversion layer (IKnutson et al.ll2010l ). 



The ultra close-in orbit of TrES-3b also makes it a great target for follow-up transit 
monitoring. It has been s uggested that plane ts with very short periods are likely to be 
falling into their host stars (I Jackson et al.ll2009l ). Some have l ooked for trarisit tim ing varia- 
tions indicative of such orbital decay (for e.g. OGLE-TR-56b lAdams et al.ll201l[ ). TrES-3b 
is similar to OGLE-TR-56b in many regards, including the orbital period (< 1.5 days), size 
(~ 1-Rjup) and impact parameter {h > 0.8). Though there has been no str ong evidence 



for T TVs for OGLE-TR-56b, it is one object that warrants long term study (lAdams et al. 



20111). Tra nsit monitoring of TrES-3b by other tea ms has so far confirmed its linear orbital 
ephemeris (ISozzetti et al.ll2009l : iGibson et al.ll2009l ). A search for additional transiting com- 
panions using the EPOXI mission has yielded no detections, since the likelihood of detecting 
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a planetary sibling in a c oplanar orbit yyith TrES-3b is lower for an outer planet given the 
inclination of the system (IBallard et al.ll201l[ ): inner planets are not expected since TrES-3b 
is already quite close to its parent star. 

In this paper we report observations of eleven transits of TrES-3b, taken as part of 
the Apache Point Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE). 
The APOSTLE program is a follow-up transit monitoring program designed to obtain high- 
precision relative photometry on known transiting syste ms in order to refine me asurements 
of system parameters and transit times (see e.g. GJ 1214b iKundurthy et al.ll201ll ). In §[2] we 
outline our observations, and in § [3] we describe the data reduction. The two sections which 
follow, § m and § O outline the transit model and the Markov Chain Monte Carlo (MCMC) 
analyzer respectively. In § [H] we present our estimates of the system parameters for TrES-3b 
and in the subsections § 16.11 and § 16.21 we present results from our study of transit depth 
variations (TDVs) and transit timing variations (TTVs). Finally, in §[7] we summarize our 
findings. 



Observations 



TrES-3 was observed by APOSTLE over a time span of 2 years between the summer of 
2009 and the fall of 2011. All obs ervations of TrES-3 we re carried out using Agile a high- 
speed frame-transfer photometer (IMukadam et al.l 120111 ) , on the ARCEI 3.5m telescope at 
Apache Point, New Mexico. The Agile CCD has no dead time, as the charge is transferred 
to an adjoining array for read-out. All observations for TrES-3 were made using Agile^s 
medium-gain, slow read-out (lOOkHZ) mode, with the charge read-out at 45sec intervals. 
The filter used for the observations was the /-band ^imita to the SDsi r filter (with 
central wavelength, Aq = 626nm, iFukugita et al.lll996l ). This observing filter is bluer than 



typical filters used for transit observations since Ag ile is a blue sensitive C CD that is affected 
by a strong fringe pattern at longer wavelengths (IMukadam et al.ll201ll ). The summary of 
the 11 r'-band observations is given in Table [H Observations were made by adjusting the 
focus on the secondary mirror to smear the stellar Point Spread Functions (PSFs) across 
multiple pixels, which minimizes the systematics caused by pixel-to-pixel wandering of the 
PSF. The long read-out (exposure time) also allowed for a greater count rate that maximized 
the signal-to-noise per image. The count rate was kept below Agile's non-linearity limit of 
~52k ADU and well below its saturation level of 61k ADU by small adjustments to the 
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telescope's secondary focus. 

The parent star of TrES-3b, GSC 03089-00929 (TrES-3) is a G type star with a Johnson 
R magnitude of 12.2. The comparison star used for the relati ve photometry was USNO-Bl.O 



1275-0332540 situated ~ 90" away, with Johnson R of 12.9 flMonet et al.ll2003f ). which was 
the next brightest star in Agile's FOV. Several studies that use detectors with larger FOVs 
than Agile typically use many (N >> 1) comparison stars. For relative aperture photometry, 
lightcurve precision is limited by the signal-to-noise achieved from aperture extraction on 
the faintest star in the set. The small FOV of Agile meant that most other comparison stars 
in the field were too faint to provide the adequate signal-to-noise on the final lightcurve. 
PSF photometry is a solution to the problem of using faint comparison stars since the stellar 
and background components can be constrained accurately. However due to the difficulty in 
modeling the complex defocused PSF of APOSTLE observations we did not opt for this route. 
We used a circular aperture for photometry (see details in § [3]) and hence the target and 
brightest companion produced the best results. Our uncalibrated differentrial photometry 
shows that TrES-3 was the brighter of the two by a factor of ~ 1.8 in the r'-band. We 
also note variability in the uncalibrated flux between TrES-3 and its reference star, with the 
maximum difference being ~20% between the highest (#7 UTD 2011-03-24) and lowest (#3 
UTD 2010-03-22) values. The observations were made over a variety of observing conditions 
(Column 'Obs. Conditions' in Tabled]). The listed nights include both complete and partial 
transits (where data were lost due to poor weather, or issues with the instrument). Some of 
these partial transits include UTD 2009-05-14 (#1), 2009-06-13 (#2) and 2010-06-02 (#10). 
Small portions of the in-eclipse data were lost for the transits on UTD 2010-03-22 (#3) 
and UTD 2011-04-27 (#8), and the night of UTD 2011-05-14 (#9) had exceptionally poor 
observing conditions (seeing > 2"). Even though these data will affect the fit, we include 
them in the analysis since they can be used to determine transit times and other system 
parameters. 



3. APOSTLE Pipeline 

The APOSTLE project used a customized data reduction pipeline, written in the In- 
teractive Data Language (IDL) to process data from Agile. The pipeline performs standard 
image processing steps like dark subtraction and fiat fielding, but also implements non- 
linearity corrections unique to Agile. The pipeline also creates an uncertainty map of the 
processed images by propagating pixel-to-pixel errors through each step of the reduction. 
In addition to the photon counting errors and read-noise from the raw images, the pipeline 
propagates the variance on the master dark and master fiat during the reduction. Errors 
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were also propagated for those pixels where the counts exceeded the non-linearity threshold 
(~52k counts) using the uncertainties in an empirically derived non-linearity correction func- 
tion. Frames where pixels inside a photometric aperture exceeded Agile^s saturation limit of 
61k were rejected. Images at the other extreme, where the stars were obscured by clouds, 
and resulted in low signal to noise measurements were also rejected {i.e. where photometric 
errors were > SOOOppm). The fraction of rejected frames per night is listed in Column (9) 
'%Rej.' in Tabled 



We used SExtractor ( IBertin &: Arnoutd Il996l ) to derive initial centroids of our de- 



focused stars. This software allows the use of a customized PSF kernel, so we used a 
'do nut '-shaped detection kernel for our defocused data. Coordinates obtained from SEx- 
tractor were then used for circular aperture photometry with the PHOT task in IRAF's 
NOAO.DIGIPHOT.APPHOT package. We derived flux estimates from a range of circular 
apertures with radii between 5-50 pixels, at intervals of 1 pixel, by simply summing the 
counts in these apertures. An outlier-rejected global median on the frame was used as the 
sky estimate, which is removed to derive the instrumental flux of the stars. To derive pho- 
tometric errors, we extracted counts from the error frames using the same centroids and 
apertures used for photometry on the target frames. The lightcurves are generated by di- 
viding the instrumental flux of the target star by the comparison star and then dividing the 
entire lightcurve by the median out-of-eclipse flux level. At this stage there may still be 
systematic trends in the lightcurve which have not been removed by reduction; for example, 
differential extinction due to airmass variation or photometric variation due to centroids 
wandering over pixels of varying sensitivities (e.g. due to small imperfections in the flat- 
fielding). Thus, for each image we also extracted a set of nuisance parameters which are used 
to compute a correction function {i.e. detrending function). For the TrES-3 data we found 
that, (i) the airmass, (ii) the global median sky, (iii) the centroid positions of the target 
star, and the sum of counts in the photometric aperture for (iv) the master dark and (v) the 
master fiat, showed trends that corresponded to trends seen in the lightcurves. The airmass 
was derived from the image headers, while the global sky and the centroid positions are 
derived from the photometry on the science frames. The sum of counts in the master dark 
and master flat are derived from photometry on the master dark and flat using the centroids 
and apertures used for photometry on the science frames. The correction function {F^or) is 
modeled as a linear sum of nuisance parameters as described by the following equation: 



k=l 



where j are the nuisance parameters, Ck are the corresponding coefficients. The index k 
counts over the number of nuisance parameters iVnus, and the index i denotes the transit 



6 



number. The detrending coefficients are chosen by minimizing the between the observed 
data (O), a model function (M) and correction function, 



here j is the index which counts over the total number of data points {N^ii), when all transits 
are stacked. During photometry we use a set of trial model parameters (based on values from 
the literature) to remove the transit lightcurve from the data. We use a linear least-squares 
minimizer to fit for the coefficients of -Fcor.i (Eq- [I]), which removes any correlated trends that 
remain in the model-subtracted lightcurve. The resulting residuals are used for selecting the 
optimal aperture from which to extract the photometry. The aperture where the scatter in 
the residuals is minimized was chosen as the optimal aperture. Smaller apertures do not 
completely sample the flux from the stars, and at larger apertures one loses signal-to-noise 
due to the fact that one accumulates fewer counts from the star, compared to sky counts. 
The optimal aperture typically fell between 5-50 pixel radii. We present the complete list of 
optimal apertures in Column (6) in Table [TJ The sizes of the typical apertures were always 
a few times larger than the typical half-width at half maximum of the stellar PSFs. The 
photometric precision, which is simply the scatter in the residuals, is also listed in Table [1] 
in Column (7). 




(2) 



Table 1: APOSTLE Observing Summary for TRES3 
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(31 
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(71 


(81 


(91 


(101 


1 


2009-05-14 


Clear 


r' 


45 


30 


777 


2% 


1.9179 


0.5618 


2 


2009-06-13 


Poor Weather 


r' 


45 


22 


610 


13% 


1.9083 


0.5135 


3 


2010-03-22 


Clear 


r' 


45 


20 


962 


< 1% 


1.7788 


0.7436 


4 


2010-05-16 


Clear 


r' 


45 


21 


877 


1% 


1.8954 


0.7769 


5 


2010-06-02 


Poor Weather 


r' 


45 


30 


703 


< 1% 


1.8670 


0.2287 


6 


2010-10-12 


Clear 


r' 


45 


24 


896 


< 1% 


1.8832 


0.7804 


7 


2011-03-24 


Clear 


r' 


45 


19 


644 


1% 


1.9809 


0.5644 


8 


2011-04-27 


Poor Weather 


r' 


45 


14 


1022 


1% 


1.8932 


0.9787 


9 


2011-05-14 


Poor Weather 


r' 


45 


15 


2960 


< 1% 


1.8887 


2.8965 


10 


2011-06-21 


Clear 


r' 


45 


26 


1221 


< 1% 


1.9227 


0.9683 


11 


2011-08-24 


Clear 


r' 


45 


19 


1240 


7% 


1.8016 


0.7597 



(1) Transit Number, (2) Universal Time Date, (3) Observing Conditions, (4) Observing Filter, (5) Exposure Time (seconds) 
(6) Optimal Aperture Radius (pixels), (7) Scatter in the residuals 

(8) % frames rejected due to saturation or other effects, (9) Flux normalization between the target and comparison star 
(10) The factor by which the photometric errors were scaled 
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Table 2: APOSTLE Lightcurve data* for TrES-3 



T# 


T-TO 


Norm. Fl. Ratio 


Err. Norm. Fl. Ratio 


Model Data 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


-0.0490181 


0.999407424 


0.000934904 


1.000000000 


1 


-0.0479764 


1.000889277 


0.000936546 


1.000000000 


1 


-0.0474555 


1.001329344 


0.000936130 


1.000000000 


1 


-0.0469347 


1.002419716 


0.000941452 


1.000000000 


1 


-0.0455226 


1.000658752 


0.000936362 


1.000000000 


1 


-0.0444809 


0.998944209 


0.000923458 


1.000000000 


1 


-0.0434392 


0.998014984 


0.000914333 


1.000000000 


1 


-0.0429184 


0.996495746 


0.000917120 


1.000000000 


1 


-0.0423975 


0.999418076 


0.000914379 


1.000000000 


1 


-0.0418767 


1.000283709 


0.000908789 


1.000000000 



*Thc data are presented in their entirety as an online-only table. 
(1) Transit Number, (2) Time Stamps - Mid Transit Times (BJD) 
(3) Normalized Flux Ratio, (4) Error on Normalized Flux Ratio 
(5) Model Data 
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The 11 transits of TrES-3b are shown in Figured] in normahzed flux ratios (with offsets 
for clarity). The plotted data result from the data reduction and model fitting processes 
described in sections § El § IH and § O respectively. The plotted data are presented as an 
online-only table. 



4. Multi Transit Quick 



We developed a transit model called MultiTransit Quick (MTQ) in 



is based on the analytic lightcurve models presented in iMandel fc Ago! 



'YTHON, which 



(120021) ■ and th e 



PYTHON implementation of some of its functions (from EXOFAST by lEastman et al.ll2012l ). 
The set of transit parameters used by MTQ include the transit duration (^t), the limb-crossing 
duration {to) and the times of mi d-transit (TO). The parameters It and to are the same as 
T and r from ICarter et al.l (120081 ). The shape of the transit is characterized by the transit 
depth (D) and the stellar limb- darkening parameters Vi and V2- The parameter D simply 
represents the maximum depth of the transit trough at conjunction, which depends on the 
ratio of the area of the disks of the planet to the star (Rp/Rl) and the ratio of the disk- 
averaged stellar intensity at conjunction I{b) to the un-obscured disk-averaged intensity at 
conjunction for an impact parameter of zero, /(O). The transit depth D is: 



D 



Km 
Ri m 



(3a) 



the t erm I(b)/I(0) can be replaced by the quadratic limb- darkening profile ( iMandel &: Ago! 
2002h to give 

^ l-«i(l 



D 



Rl 



:i-f - 



U2 ' 



(3b) 



wher e Ui and U2 are the quadratic limb-darkening coefficients (same as 71 and 72 in IMandel fc Agol 
2OO2I ) . The limb-darkening parameters Vi and V2 in MTQ are linear combinations of the 
quadratic limb- darkening coefficients ui and M2, with f 1 = mi -|- U2 and V2 = ui — U2- These 
linear combinations were used while fitting the lightcurves, since it is known that directly fit- 
ting for limb- darkening coe fficients results in s trongly anti-correlated error distributions for 
various transit parameters (IBrown et al.ll200ll ) and severely hinders the chance of Bayesian 
techniques from converging to accurate values (more in § [5]) . 

We wrote two versions of MTQ, one designed to fit transits observed with different filters 
(Multi-Filter), and the other designed to look for variations in the transit depth (Multi- 
Depth). If one makes transit observations using filters of different wavelength, the resulting 
set of transit lightcurves must be analyzed using a model that accounts for the different 
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limb-darkening profiles, and differing transit depths extant in the data. The standard set 
of parameters used for Multi-Filter version of MTQ is ^Muiti-Fiiter = {tr, tc: Dj ^p, ?^i,j...iVF5 



V2,j...Np 



Ti,,,NT}^ where Tj are the transit times and Vi and V2 are the limb-darkening pa- 



rameters described in the previous paragraph. The subscripts {...Nt and j...Np are used 
to denote multiple transits (Nt) and multiple filters (Np) respectively. For APOSTLE's 
TrES-3b data we only observed using one filter, and the number of transit was eleven. 

The second version of MTQ is designed to fit for the depths of each individual transit 
lightcurve separately. Variations in transit depth can arise for several reasons. Commonly 
invoked sources of such variations are starspots (cool photospheric regions) and f aculae (hot 



photospheric regions) (iPont et al.l l2007l : iLanza et al.l l2009l : iKnutson et al.l I2OIII . and many 
more). The appearance and disappearance of spots on the stellar surface due to rotation 
or stellar activity cycles would result in variations in the transit depth. Transits occurring 
across the spotless stellar surface will be shallower than when active regions exist on the 
visible face of the star (for cool spots). Conversely, if these active regions are hot spots, 
the transit depths would change in the opposite manner. Other sources for tr ansit depth 
variat ions include planetary o blateness, spin precession fjCarter &: W inn 2 0101) . planetary 
rings (IBarnes fc Fortneyll2004l ). and satellites (ISartoretti fc Schneideriil99a : iTusnski fc Valio 
2OIII ). The wide variety of proposed sources of transit depth variation means there may be 
several degeneracies to resolve if a transit depth variation is indeed detected. Nonetheless, 
understanding such phenomena may only be possible by establishing statistically significant 
measurements showing variable transit depths. The multi-filter capability of MTQ can be 
modified to fit for the depth of each transit as a unique parameter. In this case the set 
of parameters used is ^Multi-Depth = {tr, to, Di,„NT, ^ij-Np, V2,j...Np, Ti,„Nj.}, where A is 
now fit for each transit instead of each filter. However, the filter corresponding to each 
depth is still tracked as the model needs to convolve the limb-darkening profile to correctly 
reproduce the full lightcurve profile. For TrES-3b observations only a single filter was used. 
In addition, one transit is chosen as the "reference" transit, using which MTQ internally 
computes several transit parameters such as Rp/Ri,, a/R^, etc. We picked transit # 7 since 
it had high photometric precision, and was well-sampled. 

In the sections that follow we describe results from our attempt to fit for transit param- 
eters using both the Multi-Filter and Multi-Depth models. The following section outlines 
the Markov Chain Monte Carlo analyzer used in conjunction with MTQ. 



Transit MCMC 



Markov C] 
on exoplanets ( 


lain 


Monte Carlo (MCMC) analyzers are now the standard for modeling data 


Fore 


2OO5I: 


Holman et al. 


2006 


Collier Cameron et al. 


2007; 


Kundurthv et al. 


2011; 


Gazak et al. 


2011 


). We developed an MCMC routine called Transit MCMC (TMCMC), 



which is designed to work in conjunction with MTQ. However, it can also be used to fit 
other models and data easily. The core TMCMC routine uses the Metropolis-Hastings (M- 
H) algorithm, and is base d on i ts im plementation for astronomical data, as described by 



Tegmark et al.l (120041 ) and iFordI (120051 ). The Markov chain is computed by making jumps in 
parameter space and selecting those jumps which tend toward regions of parameter space at 
lower x^. Jumps from the j^^ step in the chain are made with the following equation: 



e 



(4) 



where d and ctq are the vectors of model parameters and their associated step-sizes respec- 
tively. The term ^(0, Ug), referred to as the proposal distribution, is a random number drawn 
from a normal distribution with a mean of and a variance of a^. It is customary to use a 
Gaussian proposal distribution but not mandatory. We modify the M-H algorithm slightly 
by allowing adaptive jump-size adjustments that optimize the sampling rate of the Markov 
Chain integrator. The desired acceptance rates are achieved by adjusting the step-size con- 
troller (/) every 100 accepted steps according to /new = W^l/ow/A'triais), where A'triais are 
the number of steps attempted for the last 100 accepted steps and is a scahng factor 
which is 225 or 434 for singl e-parameter or multi-parameter chains respectively (as noted in 
Collier Cameron et al.ll2007l ). Given sufficient time, a converged chain will have traversed 
parameter space such that the ensemble of the points accurately represent the uncertainty 
distributions of the parameters in set d. 



5.1. Markov Chains for MTQ 



For APOSTLE data sets, we explored system parameters using three different kinds of 
chains. Two of these were based on the Multi-Filter parameter set ^Muiti-Fiiter described in 
§IH First, with Fixed Limb-Darkening Coefficients, and second, with Open Limb-Darkening 
Coefficients. The term limb-darkening coefficient will henceforth be abbreviated as LDC. 
For the Fixed LDC chains (FLDC), the coefficients were sirn ply fixed to values tabulated 
for the appropriate observing filter (jClaret fc Bloemenll201ll ). For the Open LDC chains 
(OLDC), the limb-darkening parameters Vi and V2 are allowed to float. It has been shown 
that measured limb-darkeni ng coefficients frorn high precisiori stud ies are in good agree- 
ment with tabulated values ( Brown et al. 2001 ; Tingley et al.lboOGi ). Another study using 
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spectrophotometry from HST STIS showed that any variations in the estimate of transit pa- 
rameters due to inaccurate hmb-darkenin g are lower than the l-cx uncertainty, with greater 



disagreement seen at shorter wavelengths ( iKnutson et al.ll2007l ). Ground-based observations 



typically cannot achieve the level of photometric precision of space-based studies, hence 
constraining the limb-darkening from transit observations is difficult. Subtle inaccuracies 
in the fit limb-darkening profile can lead to incompati ble estimates of system parameters. 



especially at short wavelengths (IKundurthy et al.ll201ll ). In effect, the purpose of the Fixed 



LDC and Open LDC chains are simply to compare the differences in the fit limb-darkening 
coefficients to the tabulated values in the literature. The third type of Markov chain was 
run on the Multi-Depth parameter set ^Multi-Depth described in § HI APOSTLE lightcurves 
were gathered over a long time-baseline, and statistically significant depth variations seen 
in the data may help shed light on the various phenomena responsible for depth variations 
(see §11]). 

Table 3: Bounds applied for MultiTransitQuick in TMCMC 

Bounds Notes 

> Non-zero transit duration 

> Non-zero limb-crossing duration 
-D > Non-zero transit depth 

1 — 6^ > Impact parameters less than 1 (primary condition for transit) 
h/{a/Ri,) < 1 Ensures real values for orbital inclination 
< Ml < 1 Reasonable limb-darkening coefficients* 
< ^2 < 1 Reasonable limb-darkening coefficients* 
*applied when parameters were fit (OLDC chains) 



We applied bounds to several transit parameters and combinations of transit parameters 
in MTQ (as shown in Table[3]). Most of these bounds were simply to check if the parameters had 
values that were physically realistic. For most parameters it was quite rare that the MCMC 
chains got close to the bounding limits, since the chains spend most of their time near low 
regions; best-fit parameter values for the TrES-3 system are far from any physical bounds. 
When fitting for the limb- darkening however, the degeneracies proved very significant, and 
the Markov chains often strayed to unrealistic values. For example Ui and U2 took on values 
which would suggest limb-brightening rather than limb-darkening, either as a result of noisy 
data or due to gaps in the ingress or egress portions of the lightcurve. Limb-brightening 
is considered unphysical for broadband observations. Thus, after every jump (Eq. H]), the 
vector of proposal parameters Oj^i is run through a series of sub-routines in MTQ, to check 
if any of the conditions in Table [3] are violated. If any one of these conditions are violated. 
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the vector is discarded and a new vector is generated. TMCMC carries out this check until 
agreeable values emerge, and then continues with the rest of the algorithm. 



5.1.1. Executing and Analyzing Chains 

By varying the entire vector of model parameters, and applying a single step-size mod- 
ifier (Eq. H]), we run the risk of using mismatched step-sizes and hence undersampling the 
posterior distributions of some parameters. As mentioned before, well-cons tructed chains 



will p roperly sample posterior distributions given the correct acceptance rate ( iGelman et al. 
20031 ). The key to properly constructing a chain is to choose the relative starting step-sizes, 
for the parameter ensemble, such that they all roam high and low probability regions of pa- 
rameter space at roughly the same rate. Determining a reasonable set of starting step-sizes 
is done by running short exploratory Markov chains (40,000 steps) for each model parameter 
(holding all others fixed). If these e xploratory chains ha ve not stabilized to the optimal ac- 



ceptance rate of ~ 44% (as noted by lGelman et al.ll2003l for single parameter chains) at the 



end of 40,000 steps, the chain is run until this rate is achieved. The jump sizes near the end 
of stabilized exploratory chains are then used as the starting steps for the multi-parameter 
chains. These multi-parameter chains will also be referred to as 'long chains' from now on. 

For each transiting system, we ran long chains of 2 x 10^ steps from two different 
starting locations for each model scenario: Fixed LDC, Open LDC and Multi-Depth/Fixed 
LDC After completion we (1) cropped the initial stages of these chains to remove the burn-in 
phase, where the chain is far from the best-fit region, and (2) we exclude the stage where the 
chain is far from the optimal acceptance rate of 23± 5%, as noted for multi-parameter chains 



(IGelman et al.ll2003l ). We run three types of post-processing on the chains after cropping: 
(a) We compute the ranked and unranked correlations in the chains of every fit parameter 
with respect to the others. These statistics provide an estimate of the level of degeneracy 
between parameters in a given model. The next post-processing steps are two commonly 
used diagnostics to check for chain convergence, na mely (b) computing the auto-correlation 



lengths and (c) the Gelman- Rubin R-static values (IGelman fc Rubin! 119921 ). 



In a Markov chain, since a given step is only dependent on the preceding step, sections in 
the chain may show trends {i.e. are correlated). The correlation length signifies the interval at 
which sampled points in a Markov chain will be uncorrelated. The effective length is the total 
number of points in the chain divided by the correlation length. Short correlation lengths 
[i.e. large effective lengths) indicate that the MCMC ensemble represents a statistically 



signifi cant and hence more precise sampling of the posterior distribution (ITegmark et al. 



20041 ). Effective lengths > 1000 steps are commonly considered to be satisfactory. The R- 
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statistic is computed using multiple chains of the same parameter set which have different 
initial conditions. An R-statistic close to 1 (to within 10%) indicates that all chains have 
converged, cover approximately the same region of parameter space, and that the relevant 
parameter space has been sufficiently explored. Results from a chain are deemed useful if 
the auto-correlation and Gelman- Rubin conditions have been met. 



5.2. Comparing with the Transit Analysis Package 



It has been noted by ICarter fc Wind (120091 ) that transit lightcurves lacking any sig- 
nificant "defects" or artificial trends can have correlated noise buried within the overall 
scatter, which is invisible to visual examination. They test a wavelet based red-noise model 
on simulated transit lightcurves with and without artificial red-noise and find that models 
which do not fit for red-noise are subject to inaccuracies in transit parameters on the order 
of 2-3(7 and tend to have underestimated errors by up to 30%. For transit timing studies, 
poor estimates such as these are cause for concern, since smaller errors and large deviations 
from the expected time can easily lead to false claims of TTVs. The detrending routine 
within TMCMC removes long-term trends related to instrumental or other nuisance parame- 
ters. We visually examined the residuals of APOSTLE lightcurves (those observed during 
good conditions), and noted that the data may still h ave low-level corre lated noise, even 
after detrending. The Transit An a lysis Package (TAP iGazak et al.l l201ll ) implements the 
red- noise model of ICarter fc WinnI (120091 ). We ran separate fits of transit parameters using 
TAP on the detrended lightcurves from our TMCMC fit. 

The typical TAP parameter set is: 0tap = {a/R*, ^orb, {Rp/ R*)i...NF, Ti,„Nj., o-(„ute,i...NT), 
c"(red,i...Afr)}) where cr(„hite,i...AfT) ^^cl cT(red,i...AfT) ^hc white-noisc and red-noise levels for Nt 
transits, respectively. The TAP package does not fit for the period using the transit times, 
and often yields poor estimates of the period, so we fixed the period to that from TMCMC 
fits. The limb-darkening was fixed to values from the literature. The orbital eccentricity and 
argument of periastron were kept fixed at for TrES-3b. 



One must also note a pitfall of the ICarter fc WinnI (120091 ) red-noise model. Their fit- 
ting function expects the white-noise and red-noise components of a given lightcurve to 
be stationary, i.e. there are no temporal variations in the Gaussian white-noise level and 
the red-noise's power-spectrum amplitude with time. They note that more elaborate noise 
models may be required to account for the fact that real data do not conform to these re- 
quirements. Typically the scatter is increased (due to lower signal-to-noise) at the portion 
of the lightcurve where the sky brightness is higher (like data taken close to twilight). Vari- 
able observing conditions may alter not just the white-noise properties of a lightcurve but 
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also re sult in red-noise that is more complex than what can be described by ICarter &: Winn 
( 120091 ) 's wavelet model. In addition, at the 1000 ppm level, stellar variability is not under- 
stood; it is reasonable to assume that the target and comparison may each contribute to 
the correlated noise in the lightcurve. In spite of these caveats, the TAP red-noise analysis 
serves as a good secondary check to results derived using TMCMC (which does not account for 
red-noise). 



6. System Parameters 

As described in § 15.1.11 we ran three chains on lightcurves of TrES-3b: using the 
^Multi-Filter parameter set with (1) Fixed LDC, (2) Open LDC, and (3) using the ^Muiti-Dopth 
parameter set with Fixed LDC. The parameter sets are described in § |H All TrES-3b chains 
use tci tr and include the transit times Tj for the 11 transits. Since all TrES-3b observations 
were taken in the r'-band, the Fixed LDC and Multi-Depth Fixed LDC chains fit for 2 LDCs 
for the r'-band. The Multi-Depth chain was set up similarly, but with 11 free parameters 
for the transit depths. Post processing statistics and other data for these chains are listed in 
Table m The column Wfree', 'Chain Length', 'Corr. Length' and 'Eff. Length' list the num- 
ber of free parameters, the length of the cropped chain, the correlation and effective lengths, 
respectively. All chains were run for approximately 2 million steps, but about 100,000 of 
the initial steps were removed to account for "burn-in" and selection rate stabilization. The 
Open LDC and Multi-Depth chains both have quite low effective lengths indicating poor 
Markov chain statistics. Only the Fixed LDC, ^Muiti-Fiiter model satisfies the condition of a 
well sampled posterior distribution (effective length is > 1000). The final two columns list the 
goodness of fit {i.e. lowest in the MCMC ensemble) and Degrees-of-freedom (DOF) from 
the respective chain. Parameters from all chains had Gelman-Rubin Il-statistics close to 1 
indicating that the parameter space was covered evenly (though the OLDC and Multi-Depth 
FLDC chains were not sampled finely enough, based on the auto-correlation data). 



Table 4: TMCMC Chains for TrES-3 



Chain 


Model Vector 


Ntree 


Chain Length 


Corr. Length 


Eff Length 




DOF 


FLDC 


^Multi-Filter 


14 


1,900,001 


318 


5,974 


2503.67 


2575 


OLDC 


^Multi-Filter 


16 


1,900,001 


4,996 


380 


2519.14 


2573 


MDFLDC 


^Multi-Depth 


24 


1,900,001 


2,866 


662 


2692.94 


2565 



As noted previously in iKundurthy et al.l (120111 ) the r'-band limb- darkening is difficult 
to characterize using transit models. This result is seen again given the fact that the OLDC 
chain has the worst auto-correlation statistics. In addition the fact that TrES-3b has a near- 
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grazing transit makes constraining the transit shape more challenging. The Multi-Depth 
models also did not converge, which could be due to the fact that several of the transits are 
not completely sampled during the eclipse event (see Figured]). 

Table 5: TrES- 3 Parameters for ^Muiti-FUtcr 

Parameter FLDC OLDC XMt 



MTQ Parameters 



tG 


0.0210±0.0004 


0212+°-°°°'^ 


days 


tr 


0.0383±0.0002 


0.0385±0.0010 


days 




0.0251±0.0002 


0.0255±0.0003 






(0.6767) 


0.7371±0.1279 






(0.3008) 


-0.1782±0.4955 




Derived Parameters 




0.1652±0.0009 


0.1649±0.0015 




b 


0.836±0.003 


0.837±0.008 




a/R^ 


5.97±0.03 


5.911^^ 




'^orb 


81.95±0.06 


81.86l°:°« 


%deg) 


v/R^ 


28.71±0.14 


28.42i°i^ 


days~^ 




2.36±0.03 


9 9Q+0.05 
^•^^-0.06 


g/cc 


P (1.3062 days +) 


-1141±21 


-1109±22 


milli-sec 



The resulting best-fit parameter estimates are listed in Table [5] for the Multi-Filter 
models, and in Table [6] for the Multi-Depth models. These tables also list the derived 
system parameters. The tr ansformation betwe en th e MTQ parameters to the derived system 
parameters are described in Carter et al. ( 20081 ) and Kundurthy et al. (2011). Contour plots 
showing the joint probabihty distributions (JPDs) for the fit and derived parameters are 
shown in Figures [2] and [3] respectively. 

It is interesting to note that there is a clear degeneracy in the MTQ parameters as seen 
by the correlations in the posterior distributions of three parameters {D, to and tx) shown 
in Figure [21 From the impact parameter (6 = 0.84) we note, as previous studies have, that 
TrES-3b has a near- grazing transit. The correlation in t he JPDs seen in F igure [2] can be 
understood from Equation [3b] and the following equation ([Carter et al.ll2008l ): 



l-itT/tG){Rp/R. 



(5) 



Rearranging these two expressions gives an expression where ^ oc . [Carter et al.l (|2008[ ) 
have shown that when b ^ 1 [i.e. the transit is close to grazing), the fraction tc/tT rises and 
the covariances between D, tc and rapidly deviate from zero. This degeneracy may be the 
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Table 6: TrES-3 Parameters for ^Muiti-Dcpth 



Transit Depths 


Value 


Units 


Rp/R* 


Value 


Units 




0.0249±0.0007 


- 


{Rp/Ri,)i 


0.1676±0.0025 


- 


{D)2 


0.0261±0.0009 


- 


{Rp/Ri,)2 


0.1711±0.0030 


- 


(Dh 


0.0276±0.0007 




{Rpl R*)3 


0.1754±0.0024 






0.0256±0.0005 




{Rp/ i2*)4 


0.1696±0.0021 




(Dh 


0.0237±0.0010 




{Rp/Ri,)^ 


0.1638±0.0034 






0.0261±0.0006 




{Rp/Ri,)^ 


0.1710±0.0024 






0.0259±0.0005 




{Rp/R.)7 


0.1705±0.0021 




{D)s 


0.0269±0.0006 




{Rp/R*)8 


0.1734±0.0024 




{D)9 


0.0274±0.0012 




{Rp/Ri,)g 


0.1750±0.0038 




(DU 


0.0254±0.0005 




{Rp/Ri,)io 


0.1690±0.0022 




{D)u 


0.0277±0.0007 




{Rp/Ri,)ii 


0.1758±0.0025 




Other MTQ Parameters 


Parameter 


Value 


Units 


Parameter 


Value 


Units 


to 


0.0232±0.0007 


days 


tT 


0.1758±0.0025 


days 


^i(r') 


(0.6767) 






0.1758±0.0025 




Derived Parameters 


b 


0.858±0.005 




a/R^ 


5.96±0.03 




^orb 


81.72±0.09 


"{deg) 


y/R* 


28.66±0.15 


days ^ 


P* 


2.35±0.04 


g/cc 









primary reason why there is disagreement between the estimates of the system parameters 
for TrES-3b reported by various groups (see Table [7]). It may also explain the low statistical 
significance of the Open LDC and Multi-Depth/Fixed LDC chains, since Markov Chains 
take longer to converge when there are correlations between model parameters. 

System parameters agree with previously published values in the literature, as seen 
by the overlap of the uncertainties in the JPD plot (Figure [3]). The errors from TMCMC 
on the orbital inclination (zorb), impact parameter jh) a nd si^ ellar density (p*) are smaller 
than the previous best measurements of ISozzetti et al.l (120091 ) by factors of 2.5, 3.5 and 2 
respectively. However, since the TMCMC+MTQ analysis does not include rednoise anal ysis the 
errors presented in Table [6] and Figure [3] are underestimates (ICarter fc Winrul2009l ). More 
conservative constraints were placed on a subset of these system parameters using the TAP 
package. Comparisons of some parameters and their uncertainties are presented in [3 



It is clear that using TAP on the APOSTLE dataset and accounting for rednoise provides 
more conservative estimates of the system parameters. TAP errors are 10-30% smaller 
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Table 7: Comparison of Estimates of System Parameters for TrES-3b 



Parameter 


TMCMC 


TAP S09 


G09 Cll 


Units 




5.97±0.03 


5.89±0.05 5.93±0.06 


6.01±0.84 




^orb 


81.95±0.06 


81.59±0.14 81.85±0.16 


81.73±0.13 81.99±0.30 


\deg) 


b 


0.836±0.003 


0.861±0.007 0.840±0.010 


0.852±0.013 






2.36±0.03 


2.26±0.06 2.30±0.07 




g/cc 


TMCMC & TAP values are from independent analysis of APOSTLE lightcurves 

S09 - fSozzetti et al. 2009). G09 - (Gibson et al. 2009V Cll - (Christiansen et al. 2011) 



than those reported by lSozzetti et al.l (120091 ) . ISozzetti et al.l ( 120091 ) and lGibson et al.l (|2009[ ) 
account for correlated noise by scaling the photometric errors and generally have comparable 
uncertainties to the APOSTLE analysis with TAP. One must note that improved system 
parameters derived from tr ansit measurernents , like p^, caii be u sed to place constraints 
on the age of the system (ISozzetti et al.l l2009l : ISouthworthI l2010l ) , which in turn can be 
used to underst a nd th e evolutionary history of the exoplanetary atmosphere. As noted by 



Carter fc Winru (120091 ). we confirm that correlated noise in lightcurves is an obstacle for 



obtaining more precise estimates of system parameters. 



6.1. Transit Depth Analysis 

The auto-correlation data indicate that the Multi-Depth/Fixed LDC did not converge 
and hence errors for the parameters from this chain are unreliable. However, we can still 
study the results from the Multi-Depth fit assuming that the best-fit values are accurate. 
Since the Gelman-Rubin test indicate that the parameter space was fully traversed, the 
best-fit points (median values) from the chain are likely to be close to values that may have 
resulted from a converged chain. 

Figure m shows the transit depth vs. transit epoch for 11 r'-band observations of TrES- 
3b. The overall variations in the r'-band depth are 0.12% compared to the 0.02% uncertainty 
in D(^/) from the joint fit to depths (FLDC, Tabled]). The median depth value (from depth 
measurements in Table [6]) is 0.0261, compared to D(r') = 0.0251 from the FLDC chain. We 
note that several of the lightcurves are not completely sampled. We excluded lightcurve 7^ 
3, 5, 8, 9 and 11 and recomputed the median depth and scatter to be 0.0257 and 0.04% 
respectively, which is more consistent with the results from the Fixed LDC chain. However, 
since the scatter is still larger than the uncertainty in the depth measured by the Single- 
Depth Fixed LDC chain we cannot completely rule out variability. Since TMCMC does not 
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have rednoise analysis and the transit depth cannot be computed via TAP's parameter set 
we also cannot know the level at which the depth uncertainties in the Fixed LDC chains are 
underestimates. In addition, there are several facts about TrES-3 which could support the 
notion that the depth variability seen in Figure IHis due to stellar activity. Firstly, the activity 



index reported for TrES-3 classifies it as an active star (ISozzetti et al.ll2009l : iKnutson et al. 



20101). Secondly, observations in the r'-band are known to be affected by spots, due to the 
fact that the Ha line falls in this wavelength range, and spot-to-star contrast ratios are 
enhanced. However, there have been no indications of spot-crossing events in the transit 
data (a sure sign of stellar activity); this could be due to a large impact parameter and the 
greater likelihood of starspots being equatorial (assuming the stellar spin-axis is aligned with 
the sky-plane). Thirdly, the possibility that TrES-3 is variable is evidenced by the changing 
flux-ratio with its comparison star (Column 'Flux Norm.' in Tabled]); though this fact could 
just as easily be explained by variability in the comparison. 

In summary, after accounting for incomplete and poor quality lightcurves, the result- 
ing low-level variance in transit depth does not warrant a confident claim for the detection 
of variability in TrES-3. More continuous sampling of transit lightcurves, without inter- 
ruptions would provide better insight into variability. Moreover, as with other parameter 
measurements, the affect of rednoise on depth measurements needs to be studied further. 



6.2. Transit Timing Analysis 



Using Tra n sit Ti ming Variations (TTVs) t o look for additional planets was first proposed 
by lAgol et al.l (120051 ) and iHolman fc Murrayl (120051 ): in a system where only one planet 
is seen in transit, a deviation from the Keplerian period could indicate the presence of 
additional undetected planets. TTVs on the order of minutes can be pr oduced if an unseen 



companion lies close t o mean motion resonance with the transi t planet (IHolman et al.ll2010 
Lissauer et al.ll2011al : iBallard et al.ll201ll : iNesvorny et al.ll2012l ). 



We gathered t r ansit times published by ISozzetti et al.l (120091 ) , iGibson et al.l (120091 ) and 
Christiansen et al.l (120111 ). and pooled them alongside TMCMC and TAP measurements of 11 
TrES-3 transit times. The timestamps of all APOSTLE data w ere converted to BJD (TDB) 
in the customized reduction pipeline (IKundurthy et al.ll201ll ). The APOSTLE pipeline's 
time conversions have been verified by cornpariso n to the commonly used time conversion 
routines made available by lEastman et al.l (120101 ). The transit times from the literature 
were also converted to BJD (TDB) before comparisons were made. The Observed minus 
Computed (0-C) plot is shown in Figure [51 A linear ephemeris was fit to all the data using 
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the equation, 



Ti = TO + Epochs X P 



resulting in a best fit ephemeris of, 



(6) 



P = 1.306186483 ± 0.000000070 days 
TO = 2454185.9109932 ± 0.0000502 BJD 

witli a goodness of fit = 204.1 for 52 Degrees-of-freedom. Tfie large reduced clii-squared 
X^/z^ = 3.92 indicates tfiat a linear fit does not precisely fit the transit times, and there 
are signif icant timing deviat i ons in the data. The largest timing deviation of ~ 127 sec is 
from the IChristiansen et al.l ( 120111 ) data set. As previously noted, these deviations could 
either be due to underestimated timing errors or unaccounted timing systematics between 
different studies. The standard deviation of the 0-C values is ~ 43 sec, hence the largest 
timing deviation is very close to being a 3-0" outlier. Several features of the 0-C plot seem 
to indicate that these timing variation are likely due to inconsistencies b e tween the various 
methods used to derive transit times. For example, [Christiansen et al.l (l201ll )'s times are 
consistently ea rly when compared to the expected transit times [i.e. 0-Cs are all negative). 
Removing the [Christiansen et al.l (l201l[ ) transit times from the ephemeris fit results in a 
millisecond change in the fit period, but improves the goodness of fit slightly to x^/^ = 
3.5, with the scatter being 35 seconds and the largest timing deviation being 78 sec, a ~2-cr 
offset. This fact establishes that there may be discrepancies between the methods used to 
derive the transit times, and hence we cannot claim any TTVs for TrES-3. 

In order to compare the ephemerides derived with and without rednoise analysis, we 
fit for a linear ephemeris to the APOSTLE transit times from TMCMC and TAP respectively 
(presented in the bottom half of Table [8]). The difference between the periods derived for 
these subsets and the period derived from all available transit times was < 12 milli-seconds. 
The reduced x^s were 7.99 and 1.33 for the TMCMC and TAP subsets respectively, confirming 
that TAP gives more conservative errors for the transit times thanks to the red-noise analysis. 
Given the much more robust fit to a linear ephemeris from TAP's transit times we can rule 
out any TTVs in the system larger than ~ 27.3 sec (the scatter in the 0-C for transit times 
derived from TAP). 



For the case when planets are in mean-motion resonance (MMR). lAgol et al.l ( l2005l ) show 
that the analytic expression, 5tmax ~ 4^ 7+mt — )' roughly estimate the amplitude of 
the timing deviation (5tmax)- The quantities rripert, "^trans, P and j are the mass of the unseen 
perturber, the mass of the transiting planet, the orbital period of the transiting planet and 
the order of the resonance respectively. For the TrES-3 system, we can rule out possible 
system configurations given the variations seen in the 0-C result from APOSTLE's TAP 
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Table 8: APOSTLE Transit Times for TrES-3b 



Epoch 


TO (TMCMC) 


O'TO 


TO (TAP) 




(Tto 




2,400,000+ (BJD) 


(BJD) 


2,400,000+ (BJD) 




(BJD) 


597 


54965.7046437 


0.0001079 


54965.7046300 





0002000 


620 


54995.7465390 


0.0001326 


54995.7465900 





0001400 


836 


55277.8821699 


0.0002456 


55277.8822500 





0003600 


878 


55332.7425269 


0.0001033 


55332.7426300 





0002700 


891 


55349.7228375 


0.0003216 


55349.7230400 





0003800 


992 


55481.6479691 


0.0001589 


55481.6479900 





0001600 


1117 


55644.9210892 


0.0000699 


55644.9211900 





0001800 


1143 


55678.8828740 


0.0001637 


55678.8825200 





0003000 


1156 


55695.8623990 


0.0003861 


55695.8623500 





0006600 


1185 


55733.7414753 


0.0001397 


55733.7417000 





0003300 


1234 


55797.7456860 


0.0001699 


55797.7456600 





0002900 


Fit 


Period (days) 


O-p 


TO (BJD) 




(Tto 


TMCMC 


1.306186240 ± 0.000000181 


2454965.7043612 


± 0.0000770 


TAP 


1.306186489 ± 0.000000302 


2454965.7043416 


± 0.0001118 



fit. The maximum possible TTV amplitude that could be hidden within the variations seen 
in APOSTLE transit times are cttttap ~ St 

max ~ 27s ec. Usinp^ the orbital period from 
Tableland the mass of TrES-3b, mtrans = 1-91 Mjup ( iSozzetti et al.ll2009l ). we compute 
the maximum mass perturber that could exist in the TrES-3 system in the 2:1 MMR to 
be ~ O.66M0, i.e. additional planets with Mp < 0.66 M^may exist near the 2:1 MMR. At 
higher order resonances, this maximum mass (for a possible perturber) is larger. 



7. Conclusions 



I. Photometric Precision: APOSTLE monitored TrES-3b over a period of two years 
between 2009 and 2011, gathering 11 r'-band transit lightcurves. APOSTLE achieved pho- 
tometric precision between 600-1200ppm (excluding nights with poor seeing). The summary 
of observational results is presented in Table [1] 

II. TrES-3b System Parameters: From our analysis of 11 lightcurves of TrES-3b, we 
were able to confirm previous estimates of system parameters for the TrES-3 system. Our 
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estimates of derived system parameters in Table |5] show improvements (~ 10-30%) from 
previous measurements. Due to the system's near grazing transit we note that several free 
parameters from MTQ showed strong correlations in their Joint-Probabilities (see Figure 
Such correlations are not suited for rapid MCMC convergence, yet we were able to produce 
Markov chains that converged and hence are able to derive statistically significant uncer- 
tainties for syste m parameters; see results from the FLDC chain in Table [51 As previous 



attempts showed (IKundurthy et al.ll201ll ). the Markov chains where limb- darkening coeffi- 
cients were set as free parameters failed to converge, and so we report that APOSTLE did 
not achieve the photometric precision required to constrain limb-darkening coefficients for 
TrES-3. 

III. Search for Transit Depth Variations: Variations in tra nsit depth over epoc h could 



be evidence for stellar variability. TrES-3 is known to be variable (jSozzetti et al.ll2009l ). hence 
one might have expected strong variations in the transit depth with epoch. Our uncalibrated 
flux ratios between the target and comparison show significant variability, yet we cannot rule 
out variability in the comparison star as the cause. Our Multi-Depth fits also show variations 
in the transit depth over transit epoch (see Figure [Hand Table E]), however we refrain from 
making a confident assertion on the detection of stellar variability since, (a) the parameters 
from the Multi-Depth chain do not have statistically significant errors as the chain failed to 
converge, and (b) we cannot rule out inaccuracies due to the incomplete sampling of several 
transits, and (c) the near-grazing nature of this transit makes constraining the depth of 
the transit trough more challenging. In fact removing several of the incompletely sampled 
transits lowers the overall scatter in transit depth measurements. 

IV. Search for Transit Timing Variations: The transit timing precisions achieved by 
APOSTLE easily allow for the detection of TTV signals > Imin (see TTV uncertainties in 
Table [8]). This is a direct consequence of the ability of the APOSTLE program to gather 
lightcurves with high photometric precision from the ground. However, we were unable to 
detect significant timing variations for TrES-3b in our data. We find that fitting for transit 
parameters while accounting for red-noise provided more conservative error estimates on 
transit times. Transit times derived using TMCMC showed significant d eviations from a linear 



ephemeris fit. Th e same data analyzed with analyzed with TAP (jCarter fc WinnI 12009 



Gazak et al.l l201ll ) significantly reduced these deviations. The overall scatter in the 0-C 
values from our rednoise analysis was on the order of ~ 27 sec, which rules out planetary 
companions more massive than 0.66 M^near the 2:1 MMR, and larger companions near 
higher order resonances. 
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Transit times published in the hterature are derived using different techniques for fitting 
transit parameters. A proper analysis of transit times would need a simultaneous analysis of 
transit lightcurves using a transit m odel that is (1) suit ed for Bayesian inference [i.e. with 
a fairly uncorrelated parameter set, ICarter et al.l 120081 ) and (2) a transit model that can 
adequately account for red-noise in the data (like TAP, iGazak et al.l |20 111 ). The catalog of 
transit times accumulated in the literature and used in this study of TrES-3b are not derived 
from such an analysis. MCMC analysis can also be inefficiently slow given large parameter 
sets and c omplex internal checks on a g iven model. Fast MCMC routines like The MCMC 
Hammer (iForeman-Mackey et al.l |2012| ) may in the future, allow for rapid analysis using 
complex models. We provide the entire set of APOSTLE lightcurves of TrES-3h with this 
publication to allow future projects to apply improved methods to study this system. 

Extrapolations fror n Kepler planetary c andidate data seem to indicate that small plan- 
ets may be ubiquitous (IBorucki et al.ll201ll ). Interesting trends that have b een noted are. 



(1) Hot-Jupiters tend b e alone, i.e. lacking other transiting planet siblings (ILatham et al. 



201ll : ISteffen et al.ll2012l ) and (2) members of multi-plane t systems with short period planets 



(Perio d < 10 days) are more likely to be Hot-Neptunes (iLatham et al.ll201ll : iLissauer et al. 



2011bl ). This trend seems to reveal the existence of different formation pathways among 
volatile-rich planets. Hence the lack of detections TTVs in TrES-3b (a Hot- Jupiter) is con- 
sistent with Kepler^s findings. 
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APOSTLE Observations of TRES3 



1.00 



0.95 



o 

+ 0.90 



c3 



T3 

O 
N 

I 0.85 
o 



0.80 



0.75 



Ti UTD 2009-05-14 



Ta UTD 2009-06-13 



Tg UTD 2010-03-22 



T4 UTD 2010-05-16 



T5 UTD 2010-06-02 



Te UTD 2010-10-12 



T7 UTD 2011-03-24 



Ts UTD 2011-04-27 



Tg UTD 2011-05-14 



Tio UTD 2011-06-21 




Til UTD 2011-08-24 



-0.150 -0.075 0.000 0.075 

BJD - Mid Transit Time (days) 



Fig. 1. — Eleven r'-band lightcurves of TrES-3b. The vertical axis is in normalized flux 
ratio units. The horizontal axis shows time from the mid-transit time in days, computed by 
subtracting the appropriate mid-transit time for each transit from the best-fit values in the 
Fixed LDC chain. 



- 28 - 



JPD - TrES-3 FLDC 




0.0206 0.0219 0.0378 0.0385 2.00e+04 
tc tx Number 



Fig. 2. — Plots of the Joint Probability Distributions (JPD) of parameters from the Fixed 
LDC chains, showing that due to the system's near-grazing transit the parameters chosen in 
^Multi-Filter show Correlations, unlike the cases for other systems discussed in this work. Table 
E] gives units. 
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JPD - TrES-3 System Parameters 




Fig. 3. — Plots of the Joint Probability Distributions (JPD) of derived system parameters 
from the Fixed LDC chains. Parameter estimates available in the literature are overplotted. 
Table [5] gives units. 
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Fig. 4. — The transit depth D as a function of transit number for r'-band observations 
of TrES-3b. The sohd horizontal and dashed hnes represent the best-fit value and errors 
respectively for D from the Fixed LDC TMCMC fit. The dotted line is the weighted mean of 
transit depth values from the Multi-Depth Fixed LDC chains. 
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Fig. 5. — The Observed minus Computed Transit Times for TrES-3b. Values from APOS- 
TLE's TMCMC fit, TAP and the hterature are plotted. The horizontal axis represents the 
transit Epoch. The zero-line ephemeris is described in § 16.21 



